extern doublereal diff_(doublereal *x, doublereal *y);
